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Background: Atelosteogenesis type I (AO-I) is a rare lethal skeletal dysplastic disorder 
characterized by severe short-limbed dwarfism and dislocated hips, knees, and elbows. 
AO-I is caused by mutations in the filamin B {FLNB) gene; however, several other genes 
can cause AO-like lethal skeletal dysplasias. 

Methods: In order to screen all possible genes associated with AO-like lethal skeletal dys- 
plasias simultaneously, we performed whole-exome sequencing in a female newborn hav- 
ing clinical features of AO-I. 

Results: Exome sequencing identified a novel missense variant (c.517G>A; p.Alal73Thr) 
in exon 2 of the FLNB gene in the patient. Sanger sequencing validated this variant, and 
genetic analysis of the patient's parents suggested a de novo occurrence of the variant. 

Conclusions: This study shows that exome sequencing can be a useful tool for the identi- 
fication of causative mutations in lethal skeletal dysplasia patients. 
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INTRODUCTION 

Lethal skeletal dysplasias commonly include thanatophoric dys- 
plasia, achondrogenesis, osteogenesis imperfecta type 2, con- 
genital hypophosphatasia, and atelosteogenesis (AO), and can 
be distinguished by four key features: bone mineralization, frac- 
tures, macrocranium, and trunk length [1]. AO is a rare lethal 
skeletal dysplasia. It represents a heterogeneous group of disor- 
ders with the following features: incomplete ossification of long 
bones and markedly short, hypoplastic distal femora and hu- 
meri resulting in rhizomelic dwarfism [2]. Rare types of lethal 
skeletal dysplasias, including AO, are difficult to diagnose accu- 



rately owing to their clinical heterogeneity [3, 4]. 

AO is subdivided into types I (AO-I), II (AO-II), and III (AO-Ill). 
AO-I and AO-Ill are associated with FLNB mutations [5], while 
AO-II is caused by SLC26A2 {DTDST) mutations [6]. The FLNB- 
related disorders constitute a spectrum of phenotypes ranging 
from spondylocarpotarsal synostosis syndrome and Larsen syn- 
drome at the mild end, to AO-I, AO-Ill, and boomerang dyspla- 
sia at the severe end [7, 8]. Besides FLNB, many genes, includ- 
ing /\/_P/_, ARSF^ COLlAl/2, C0L2A1, COLllAl/2, CRJAP, 
DYNC2HI FBP, FVC, FVC2, FGFR3, IFT80, IFT122, LFPRFl, 
NFKl PFX7, PPIB, SLC26A2, SLC35D1, SOX, TRlPll, and 
y\/DR35, have been identified to be involved in AO-like lethal 
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skeletal dysplasias [9]. 

Several recent studies have described how exome sequenc- 
ing has allowed the identification of causative mutations in many 
heterogeneous and complicated disorders [10, 11]. Confirming 
a definitive molecular diagnosis is often difficult and costly in 
skeletal dysplasia because patients may be in a critical condi- 
tion, in addition to having a phenotypically complex presenta- 
tion. This clinical and genetic diversity reflects the number of 
causative genes, and molecular diagnosis of such disorders by 
conventional Sanger sequencing, usually an expensive and 
time-consuming gene-by-gene screening process, is sometimes 
impractical [10]. Exome sequencing allows simultaneous analy- 
sis of all coding sequences of genes. The advantage of exome 
sequencing is that it does not require a prior knowledge of the 
genes responsible for a disorder. In this study, we performed 
exome sequencing of a newborn patient with clinical features of 
AO-I and successfully identified a de novo FLNB mutation. 

METHODS 

1. Patient 

A female baby was born as the first of twins of non-consanguin- 
eous healthy parents at 27 weeks of gestation via emergent ce- 
sarean section due to fetal distress. Both her parents were 31 yr 
old, and her mother's obstetric history was gravida 1, para 0. 
She was impregnated with twins through in vitro fertilization due 
to infertility. Prenatal ultrasonography at 25 weeks of gestation 
revealed that one fetus had polyhydramnios and multiple skele- 
tal anomalies consistent with skeletal dysplasia, such as very 
poorly ossified femora and humeri with severe shortening of all 
extremities. 

The patient had respiratory difficulty from birth and required 
mechanical ventilation support. It was very difficult to advance 
the endotracheal tube into the trachea after passing through the 
larynx. A pulmonary surfactant was instilled to treat the respira- 
tory distress syndrome of prematurity. Oxygenation and ventila- 
tion were not effective even though high-frequency oscillatory 
ventilation was applied. She died of respiratory failure at three 
hours after birth, possibly secondary of laryngeal stenosis and 
pulmonary hypoplasia. 

The clinical features of the patient indicated AO. However, 
other AO-like skeletal dysplasias could not be excluded. As con- 
ventional gene-by-gene sequencing is too costly and time-con- 
suming, exome sequencing, allowing simultaneous analysis of 
multiple genes, was performed as a post-mortem genetic study 
with the written informed consent of the parents. 



2. Whole-exome sequencing 

Genomic DNA was extracted from peripheral blood leukocytes 
by using a Wizard Genomic DNA purification kit (Promega, Mad- 
ison, Wl, USA), according to the manufacturer's instructions. 
Exome sequencing was performed by DNA Link (DNA Link Inc., 
Seoul, Korea), and the bioinformatics services were provided by 
Samsung SDS (Samsung SDS, Ltd., Seoul, Korea). The SureSe- 
lect Human All Exon 50 Mb kit (Agilent Technologies, Santa 
Clara, OA, USA) was used for in-solution enrichment of coding 
exons and flanking intronic sequences, following the manufac- 
turer's standard protocol. Adapted sequences for the lllumina 
HiSeq2000 sequencing system (lllumina Inc., San Diego, OA, 
USA) were ligated, and the enriched DNA samples were sub- 
jected to standard sample preparation for the HiSeq2000 instru- 
ment The Burrows-Wheeler alignment (BWA) was used to align 
sequence reads to the human reference genome (hgl9), and 
variants were called using the GATK software package [12-14]. 
The single nucleotide polymorphism (SNP) and short indel can- 
didates were identified at nucleotide resolution. These variants 
were annotated by ANNOVAR (version 2011 Jun 18) [15] to fil- 
ter SNPs reported in the dbSNP database (build 135) [16] and 
the 1000 Genomes Project (http://1000genomes.org). 

Polyphen-2 (http://genetics.bwh.harvard.edu/pph2/) and 
SIFT (Sorting Intolerant From Tolerant) (http://sift.jcvi.org/) were 
used to predict the consequences of the missense variants on 
protein function. The candidate variant identified by exome se- 
quencing was confirmed by using conventional sequencing. 

RESULTS 

1. Patient 

The patient's weight was 685 g (<10th percentile), height was 
22 cm (<10th percentile), and head circumference was 24 cm 
(10-25th percentile) at birth. This is in contrast to the second in- 
fant twin, who had a weight of 925 g (25-50th percentile), height 
of 36 cm (50th percentile), and head circumference of 26 cm 
(50-75th percentile). The affected infant also had a depressed 
nasal bridge, hypertelorism, micrognathia, and low-set ears at 
birth. Her thoracic cage was small, and the abdomen was protu- 
berant. She had markedly short limbs, with talipes equinovarus 
deformities and spatulated short fingers. Radiographs showed 
incomplete ossification and hypoplasia of the vertebrae, humeri, 
femora, tarsals, phalanges, and pelvis. Radiographs also re- 
vealed multiple skeletal abnormalities. The clavicles were rela- 
tively elongated, the thorax was small and bell-shaped, and sev- 
eral vertebrae had scoliosis and coronal clefts. The humeri were 
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Fig. 1. Radiograpli sliowing incomplete ossification and hypoplasia 
of the vertebrae, humeri, femora, tarsals, phalanges, and pelvis. The 
clavicles are relatively elongated (A), the thorax is small and bell- 
shaped (B), and the vertebrae have scoliosis and coronal clefts (C). 
The humeri are severely shortened and distally hypoplastic (D), the 
fibulae are completely absent, and both ankle joints, knee joints, 
and elbow joints are dislocated (E). 



All variants 
39,772 



1. Exclusion of known SNPs 



2,851 



2,407 



725 



> 

^ 2. Exclusion of variants with MAF> 1% in 1,000 Genomes Project 

^ 3. Exclusion of non-genic intronic and synonymous exonic variants 

^ 4. Screening of 23 genes through literature survey identified FLNB as candidate gene 



Fig. 2. Data analysis algorithm applied to the exome sequencing 
data. A data analysis algorithm was used to filter all the single nu- 
cleotide polymorphisms (SNPs) detected using exome sequencing, 
with numerous variants excluded at every filtering step. 
Abbreviation: MAF, minor allele frequency 

severely shortened and distally hypoplastic. The ulnae, radii, and 
tibiae were bowed. Both fibulae were completely absent. The 
ankles, knees, and elbows were dislocated (Fig. 1). 

1. Data analysis 

A mean coverage of 47.1 was achieved, and 86.6% of the tar- 
geted bases were read >10 times by exome capture and se- 
quencing. Of the 39,772 identified SNPs, the pathogenic variant 
was prioritized by using the following steps: After the exclusion of 
known dbSNP variants and variants with minor allele frequency 
(MAF) >1% in the 1000 Genomes Project, 2,407 variants re- 
mained, from which the non-genic intronic variants and synony- 
mous exonic variants were excluded. For the remaining 725 



chr3:58062997 FLBN NM_001457.3 exon 2 c.517G> A (p.Alal73Thr) 

Fig. 3. Integrative Genomics Viewer snapshot of the novel FLNB 
mutation (NM_001457.3: c.517G>A; p.Alal73Thr) identified by 
exome sequencing. 



variants, 23 genes identified through literature survey, including 
PubMed, UpToDate, and GeneReviews, were screened for de- 
termining the relevance of the remained vahants (Fig. 2). Conse- 
quently, we found one non-synonymous vahant (NM_001457.3: 
c.517G>A; p.Alal73Thr) in the FLNB gene (Fig. 3) that was 
also not reported in the Human Gene Mutation Database. This 
variant was predicted as 'probably damaging' by Polyphen-2, 
with a score of 0.998, and as 'not tolerated' by SIFT, with a prob- 
ability of 0.01, and was regarded as the best candidate for fur- 
ther evaluation. 

3. Validation of the variant 

Sanger sequencing of the FLNB exon 2 of the patient, both her 
parents, and her twin sister was performed, which revealed that 
the c.517G>A (p.Alal73Thr) variant was present only in the 
patient, indicating that the mutation occurred de novo. When a 
control group of 100 alleles from 50 healthy subjects was exam- 
ined by using exon 2-targeted sequence analysis, no carriers 
were found with the c.517G>A variant. 

DISCUSSION 

Mammals have three filamin genes: FLNA, FLNB, and FLNC. 
Filamin (Fin) A, B, and 0 comphse a family of three actin-bind- 
ing proteins, each shahng the common features of an N-termi- 
nal actin-binding domain followed by an immunoglobulin-like 
repeats domain that contains the receptor-binding region at the 
C-terminus [17]. The proteins are highly homologous, can inter- 
act to form both homodimers and heterodimers, and likely share 
similar functions across vahous organ systems. FInA is predomi- 
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nantly found in the brain and blood vessels, FInB in bones, and 
FInC in muscles. Insight into the function served by these pro- 
teins, however, has been limited to observations based on hu- 
man diseases, their expression patterns, and isolated cellular 
and molecular culture studies, largely owing to the fact that null 
FLNA and FLNC mice are embryonic lethal [18, 19]. 

FLNB mutation causes a range of skeletal disorders, from mild 
forms such as spondylocarpotarsal syndrome and Larsen syn- 
drome, to severe forms such as AO-I, AO-Ill, and boomerang 
dysplasia [7]. The most significant determinant of lethality is the 
degree of pulmonary hypoplasia. Atelosteogenesis refers to 'in- 
complete ossification,' and is characterized by defective endo- 
chondral ossification resulting in rhizomelic dwarfism [20]. AO-I 
is lethal during fetal life or early neonatal period [5]. Boomerang 
dysplasia is lethal and shows more severe defects in ossification 
than AO, and is characterized by the complete absence of ossifi- 
cation in certain limbs and vertebrae, and hypoplastic boomer- 
ang-shaped femora [4, 21]. Pulmonary hypoplasia, severe laryn- 
geal narrowing, severe tracheal narrowing, or tracheomalacia 
cause certain kinds of respiratory failure and early neonatal death 
[22]. Luewan et al. [23] reported that the combination of pulmo- 
nary hypoplasia and tracheobronchomalacia is a major cause of 
death in F/_A/B-related disorders. Larsen syndrome, AO-I, AO-Ill, 
and boomerang dysplasia are inherited in an autosomal domi- 
nant manner, and spondylocarpotarsal syndrome is inherited in 
an autosomal recessive manner [7]. 

As in the present case, most lethal F/_A/B-related disorders 
are caused by de novo mutations; thus, the risk of recurrence in 
subsequent pregnancies is very rare. The majority of the re- 
ported mutations are in exons 2-5, and mutations associated 
with the most severe phenotypes almost invariably occur in ex- 
ons 2-5 [21]. The actin-binding domain of FInB consists of two 
calponin homology domains (CHI and CH2) in its N-terminus 
[21]. Most of the causative mutations of AO-I are located in ex- 
ons 2 and 3, which encode CH2 [21]. The novel c.517G>A 
(p.Alal73Thr) variant in the present case occurs in exon 2 and 
is predicted to interrupt actin binding. 

The c.517G>A (p.Alal73Thr) variant has not been reported 
in the Human Gene Mutation Database; however, a G to T sub- 
stitution at 518th nucleotide, resulting in an alanine to valine 
substitution at the 173th amino acid, has been registered as a 
causal mutation of atelosteogenesis. A patient with the p.Ala- 
173Val mutation was reported to have vertebral abnormalities 
and carpal and pharyngeal abnormalities. Mutations in FLNB 
produce diverse phenotypes depending on their nature and lo- 
cation [24]. A different nucleotide substitution within a codon 



known to harbor a causative mutation, such as the one reported 
here, is likely to affect the phenotype in a similar manner. There- 
fore, the location of the present variant, the control study, and in 
silico analyses suggest that the novel variant in the present case 
is a pathologic mutation of AO-I. 

The diagnosis of F/_A/B-related disorders is based on clinical 
features, which might not always be conclusive, especially for 
patients in critical condition. Since exome sequencing allows 
the simultaneous analysis of all coding sequences of genes 
without a priori knowledge of their contribution to diseases, it 
can be a useful, time-saving, and cost-effective tool for provid- 
ing extended sequence information, not only for the previously 
known disease-related genes, but also to identify novel disease- 
associated genes. 

In conclusion, we presented a patient with AO-I who had a 
novel de novo heterozygous missense FLNB mutation (c.517G> 
A; p.Alal73Thr). We also showed that exome sequencing analy- 
sis can be a useful tool to identify causative mutations in patients 
with lethal skeletal dysplasias. 
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